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Abstract 



' For each n in N, let X n = [{X n )jk\^ k=1 be a random Hermitian ma- 

trix such that the n 2 random variables y/n(X n )a, v / 2ni?e((X n )j J )j < j, 
\^2nIm((X n )ij)i < j are independent identically distributed with com- 

(1) (r) 

mon distribution \x on R. Let X„ , . . . , X„ be r independent copies of 
X n and (sci, . . . x r ) be a semicircular system in a C*-probability space, 
^-j. ' Assuming that \x is symmetric and satisfies a Poincare inequality, we 

show that, almost everywhere, for any non commutative polynomial 
p in r variables, 

lim \\p(X^,...,X^)\\ = \\p(x 1 ,...x r )\\. (0.1 

, n > +oo 

J ; We follow the method of [5] and ^5] which gave ()U.1|) in the Gaussian 

(complex, real or symplectic) case. We also get that (|().1|> remains true 
£ . when the An ^ are Wishart matrices while the Marchenko-Pastur 

J> ! distributed. 
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1 Introduction 



In the 90 's, Voiculescu jTHj introduced a random matrix model for a free semi- 
circular system. He showed that if we take r independent random matrices 
(Xn )i=x,...r, distributed as GUE(n, ~), then, they are asymptotically free, 
that is, for every non commutative polynomial p in r variables, 

E[tr n p(X«,...,xW)] — r(p(x 1 ,...x r )) (1.1) 

n ► oo 

where tr n stands for the normalized trace on M n (C) and ( free 
family of semicircular variables in some non commutative probability space 
(£>, r). The result (jl.lj) holds true for a family of iid Wigner matrices and is 
proved by Dykema in 

In a recent paper, Haagerup and Thorbj0rnsen [§| proved a strong version of 
(II. lj) . in the GUE case, namely a convergence for the operator norm: 

lim \\p(xU,...,XP)\\ = \\p(x 1 ,...x r )\\a.s. (1.2) 

n — > +oo 

which led to the proof that Ext{C* ed (F 2 )) is not a group. 
Schultz [13] obtained the same result for Gaussian random matrices in the 
real case (GOE) and in the simplectic case (GSE). Our aim is to extend (jl.2j) 
in the case of an independent family of Wigner matrices on one hand and in 
the case of Wishart matrices on the other hand . Note that the special case 
r = 1 gives the well known convergence of the largest eigenvalue of to 
the right boundary of the support of X\ (see [3J for the Wigner case and [7] 
for the Wishart case; see also j2] and the references therein). 

Our approach is very similar to that of jU] and Therefore, we will 
recall the main lines of their proofs. First, in proving ()1.2)) . the minoration 

lim inf | \ P (X£\ . . . , X^) 1 1 > | \p(x h ...x r ) \ \ a.s. 

n > +oo 

comes rather easily from an a.s. version of (jl.l)) (obtained in for the 
GUE case and proved in Section 6 of [T^] for the GOE case) (see Lemma 7.2 
in [S]). So, the main difficulty is the proof of the reverse inequality: 

lim sup | b(X«, . . . , XW) 1 1 < | |p(xi, ...x r )\\ a.s. (1.3) 

n > +oo 

In the following, we sketch the main steps in the proof of (jl.3j) . 

Step 1: A linearisation trick (see [9J, Section 2 and Proposition 7.3) 
In order to prove ()1.3)) . it is sufficient to prove: 
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Lemma 1.1 For allm G N, all self-adjoint matrices do, . . . , a r x of sizemxm 
and all e > 0, 

r r 

Sp(a <g> l n + ^ a i ® X n\ U )) C S P( a ® lfi + ^ fli ® ~ C ' 6 t ( 1 - 4 ) 

i=X i=l 

eventually, as n — > oo a.e. in uj . Here, sp(T) denotes the spectrum of the 
operator T and l n the identity matrix. 

The analysis of the spectrum of S n := a <g> l n + X][=i a « ® -^n is done, using 
the Stieljes transform 

G n (X) = E[(id m <g>tr n )[(A<g>l ri -S' n )~ 1 ]], A G M m (C), im(A) positive definite. 

(1.5) 

The proof of (jl.4|) requires sharp estimates of the rate of convergence of 
Grn(A) to G?(A) := (id m <g> r)[(A <g> 1 B - s)' 1 } (of order 1/n 2 ) where s = 

Step 2: In the GUE case, Haagerup and Thorbj0rnsen [5] obtains the fol- 
lowing estimate 

||G„(A)-G(A)||<^. (1.6) 
In the GOE, GSE cases, Schultz [13] gets an extra term of order 1/n, namely 



|G„(A)-G(A)-^||<^ (1.7) 



for some functional L. 



Step 3 From the previous step, it is shown in section 6 of [0] that 

E[(ti m ®tr n ){tp{S n ))] = (tr m ®r)(<p(s)) + 0{\) (1.8) 



n 2 ' 



for ip smooth with compact support, and 



E[(tr m ®tr n )(v?(Sn))] = 0(4) (1-9) 

for smooth, constant outside a compact set and such that supp((p)nsp(s) = 
0. In the GOE case (resp. GSE case), Schultz proved in section 5 of that 

E[(tr m ®tr n ){tp{S n ))] = (tr m ®r)(^(s)) + -A(cp) + 0{\) (1.10) 



1 By a density argument, we can also assume that the matrices ai are invertible. 
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where A is a distribution with compact support in sp(s) with Stieljes trans- 
form 

/(A)=tr m (L(Al m )),AeC\R. 
Therefore, ()1.9|) still holds for <p with supp(f) fl sp(s) = 0. 

Step 4 (jl.9j) . combining with a Gaussian variance estimate, yields (by a 
standard application of the Borel Cantelli lemma), 

(tr m ®tr n )l F (S n ) =0(n~ 4/3 ) 

for F = {t E E, sp(s)) > e} which leads to (Oj) . 

The main difficulties in the generalization of the above to Wigner or 
Wishart matrices arise in step 2. Indeed, we don't have the gaussian integra- 
tion by parts' formula anymore. Our approach is inspired by the work of [T2] 
where they use a Taylor expansion (see Lemma l4.1|) extending the gausian 
integration by parts' formula. The remainder of the proof can be completed 
essentially as in the GOE/GSE case. Hence, in this paper, we shall focus on 
the obtention of such a master inequality 

\\G n (X) - G(X) - -L(X)\\ = 0(- 2 ) 

in the case of a family of Hermitian matrices with symmetric iid entries 
satisfying a Poincare inequality, as well as in the case of Wishart matrices; 
we just give some hints when the computations are similar to that of jU], [To] . 

The paper is organized as follows. In section 2, we introduce notations 
and preliminaries which will be of basic use later on. In section 3, we describe 
the proof of (jl.8|) and (jl.lUj) proved respectively in and [J3] in order to 
make clear the validity of the method in our general framework we state in 
section 4 (for the Wigner case) and section 5 (for the Wishart case). 

2 Notations and preliminaries 

This section may contain some definitions already used in the introduction 
but we choose to gather all the notations in this section for the reader's con- 
venience. To begin with, we introduce some notations on the set of matrices. 
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M p (C) is the set of p x p matrices with complex entries, M p (C) sa the 
subset of self-adjoint elements of M P (C) and l p the identity matrix. In 
the following, we shall consider two sets of matrices with p — m (m 
fixed) and p = n with n — > oo. 

Tr p denotes the trace and tr p = ^ Tr p the normalized trace on M p (C). 

||.|| denotes the operator norm on M p (C) and ||M|| 2 = (Tr p (M* M)) 1 / 2 
the Hilbert-Schmidt norm. 

Let (Eij)2j =1 be the canonical basis of M n (C) and define a basis of the 
real vector space of the self-adjoint matrices M n (£) sa by: 

'.a = i: jr ^<3<n 

e jk = —j={Ej k + E kj ), 1 < j < k < n 

fjk = ~^ E ik - E kj), 1 < j < k < n 



- For a matrix M in M m (C) <S> M n (C), we denote by 

Mij := {id m ® Tr n )(M(l m ® E 3l )) G M m (C), l<i,j <n 

and 

a>/3 fcf := (Tr m ^rf n )(M(E /3iQ ® l n )) G M„(C), 1 < a,/3 < m 

where (E a ^) is the canonical basis of M m (C). 
We now define our matrix model and the random variables of interest. 

- {Xn \ • • • , Xn^)i=i,... r is a set of iid random matrices in M n (C) sa , whose 
distribution will be specified later (matrices in GUE or GOE in section 
3, Wigner matrices in Section 4, Wishart matrices in section 5). 

- For a given family a , . . .a r in M m (C) sa , we define the random variable 
S n with values in M m (C) <g> M„(C) by: 

r 

P =i 
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and s e M m (C) <g> B by 

r 

s = a <8> Ib + ^ a p (g) a; p (2.2) 
P =i 

where the (xi)i=i,... r is a free family of self-adjoint operators in a C* 
probability space (£>, r) with a faithful state r, whose distribution will 
be specified in the different cases (semi-circular in sections 3 and 4 or 
distributed as the Marchenko-Pastur distribution in section 5). 

- For any matrix A in O where 

O :={\e M m (C)\Im(X) is positive definite}, 

we define the M m (C) valued rv: 

H n (X) = (id m <g) tr n )[(A ® l n - S n )-\ (2.3) 

G n (X) = E[H n (X)} (2.4) 

and 

G(X) = {id m ® r) [(A <g> 1 B - s)" 1 ]. (2.5) 
For A G C \ t, we also define 

# n (A) = tr ro (G„(Al m )) 

and 

51(A) = tr m (G(Al m )). 

We end this preliminary by recalling some properties of G(X) and of the 
resolvent (A (g> l n — S^) -1 of the matrix S n . First, one can easily see that for 
any A and A in M m (C) such that Im{X) and Im(X ) are positive definite, 

(A®l B -s)~ 1 -(A'®l B -s)- 1 = (A®1 b -s)~ 1 (A'-A)(A'®1 b -s)~ 1 . (2.6) 



Lemma 2.1 Let A in M m (C) s^c/i £/jat Im(X) is positive definite. Then 

|| (A (8) ls-s) _1 || < IIMA) - ^ and ||C(A)|] < ||Jm(A) _1 ||. (2.7) 
Moreover, G(X) is invertible and 

HG^A) -1 ]! < (||A|| + ||s||) 2 ||im(A) _1 ||. (2.8) 
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We refer the reader to section 5 of [H] for a proof of (|2.8|) . 

Lemma 2.2 Let A in M m (C) such that Jm(A) is positive definite, then 

\\(\®i n -s n )- l \\ < HMA)- 1 !!, (2.9) 

Vl<k,l< n,\\(\ <g> 1„ - S n )^\\ < ||Jm(A) _1 ||, (2.10) 

and for p > 2, 

1 n 

- IK*® ^ - 5 n)wT < ^IIMA)- 1 !^ (2.11) 

where, in the first inequality, \\.\\ denotes the operator norm in M m (C) (g) 
M n (C) C m M m (C) in t/ie others) and C m a constant depending only on m. 

For a Hermitian matrix M, the derivative w.r.t M of the resolvent R(z) = 
(z — M) _1 satisfies: 

R' M (z).A = R(z)AR(z) for all Hermitian matrix A. (2-12) 

Sketch of Proof: We just mention the proof of (|2.11j) . From (|2.1()|) . it's 
enough to consider the case p = 2. 

Let us denote G (n) = (A ® l n - 5 n )" 1 G M m (C) ® M n (C). Since the operator 
norm is smaller than the Hilbert-Schmidt norm, 

n 1 n m 

- 2^ W^ki II < - 2^ l^ G w I 
n 

< m||G (n) (G (n) )*|| < m\\Im(\)~ l \\ 2 . 
where the last inequality follows from ()2.9|) . □ 

In the sequel, we shall denote by any polynomial of degree k whose coef- 
ficients are positive and by C or K any constant; P^, C or K can depend on 
the ai, I — 1, . . . , r, and may vary from line to line. 
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3 Main ideas in the proofs of ( Qj ) and ( TTTTUft 
from PJ| and | |15lj 

3.1 Estimate of \\G n - G\\ in [9j 

Let us recall the main ideas of in the estimation of ||G n (A) — G(A)||. In 
lemma 5.4 of 0, Haagerup and Thorbj0rnsen observe in one hand that the 
matrix-valued Stieljes transform of s satisfies, for any A in O, 

r 

J2 a i G W a i + (ao-A) + G(A)" 1 = 0. (3.1) 
i=i 

In the other hand, using the Gaussian integration by parts formula, they 
establish the analogue of (j3.1J) satisfied by H n (X) ("Master equation", Lemma 

2 0): 



E 



0. (3.2) 



J2 a P H nWa P H n (X) + (oq - X)H n (X) + l r 
_p=i 

Then, using the Gaussian Poincare inequality to get an estimate of the vari- 
ance of H n (X), they deduce from ()3.2|) the "Master inequality" (Lemma 3 in 

0): 

|| Va i G n (A)a,G n (A) + (a - A)G(A) + l m || < ^|| MA)- 1 )) 4 . (3.3) 

Moreover, the authors prove that G n (X) is invertible for any A in O and they 
give an upper bound of the norm of its inverse (Proposition 5.2 9 ) 

HG^A)- 1 !! < (IIAH + JOIMA)- 1 !!. 

Hence, they deduce from (13. 3|) that, for any A in O, 

r 

\\a + J2 a i G nW^ + G n (A)- 1 -A|| < ^(H/mCA)- 1 !!, ||A||), (3.4) 
i=i 

where 



/ n (||/m(Ani, ||A||) = ^(||A|| +^) 2 ||/m(A)^|| 5 . 



n 
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Further, they set 

r 

A n (A) = a + ^a i G n (A)a i + G n (A) _1 
i=\ 

for any A in O. ()3.4|) can be rewritten 

||A n (A)-A||</„(||/m(A)-l,||A||). (3.5) 

The authors define 

0\ n = {\eOJ n {\\Im{\y%\\\\\)<^} 

where 

e(A) := 1 



\\Im{\)- x \\' 
(13. 5j) implies that, for any A in O n: 

ImK[x) £ 2p^rF! lm (3 ' 6) 

and that in particular A n (A) belongs to O (see Lemma 5.5 j3j). Consequently, 
applying (|3.1|) . they get that, for any A in O n , 

r r 

a + J2^G(A n (\)) ai + ^(A^A))- 1 = a + Y,a t G n (\)a t + G n (A) _1 . (3.7) 

1=1 8=1 

In proof of (b) Proposition 5.6, Haagerup and Thorbj0rnsen show that (j3.7J) 
implies that 

G n (X) = G(A n (A)) (3.8) 

for any A in 0" n : = {A G C^, e(A) > yj2 YJ i=1 INI 2 }- Using that i ^ 
/ n (t _1 , t)t _1 is a continuous strictly decreasing function from ]0;+oo[ onto 
]0;+oo[, they show in proof of (a) Proposition 5.6 9J that O n is an open 
connected subset of M m (C). Thus, by the principle of uniqueness of analytic 
continuation, ()3.8|) still holds for any A in O n . Thus, for any A in O n , they 
get that 

||G n (A)-G(A)|| < ||G(A n (A))-G(A)|| 

< ||/m(A n (A))- 1 ||||A-A n (A)||||/m(A)- 1 || 

< 2/ n (||/m(A)- 1 ||,||A||)||/m(A)- 1 || 2 
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where the last inequality comes from ([3.5)1 . (|3.6|) . Now, if A belongs to 0\0' r 
they note that 



||G n (A)-G(A)|| < 2||/m(A)- 1 || 

< ^(IIM^^iuiaidiima)- 1 !! 2 

since 

^ < / n (||/rn(A) _1 ||, ||A||)||Jm(A) _1 ||. 
Finally, for any A in O, 

\\G n (X) - G(X)\\ < 4/ ri (||/m(A)- 1 ||,||A||)||/m(A)- 1 || 2 

= ^(IIAH+^IIMA)- 1 !! 7 . (3.9) 

3.2 Estimate of \\G n -G- ±L|| in jig 

In the GOE term of order 1/n appears in the Master equation so that 

the estimate of HG^A) — G(A)|| Schultz makes by sticking to the previous 
proof of is of order 1/n. Nevertheless, a further study (we will describe 
in our general framework in section 4) gives her the sharper estimate 

||G n (A) - G(X) - -L(X)\\ < 1(||A|| + K) s P 13 (\\Im(\y l \\) (3.10) 
n n 

for any A such that ImX positive definite or negative definite. 



3.3 From Step 2 to Step 3 

^From the previous estimates ()3.9|) and ()3.10|) . Haagerup, Thorbj0rnsen and 
Schultz immediately get that, for any A in C \ M, 

l^(A)| < ^(||A]| +ir)"^(|7rr i (A)- 1 |) (3.11) 

where 

-in the GUE case 

r n (A) = g n (X) - g(\) , a = 2 k = 7. 
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-in the GOE case [S] 

r»(A) = g n (X) - #(A) - -tr m (L(Al m )) , a = 8 fc = 13. 

n 

Since and s are selfadjoint, by the spectral theory, there exist unique 
probability measures \i n and \i on R such that 

= E[(tr m ®tr B )(y>(5 ft ))] 



<^d/i = (tr m (g>r)(v?(s)). 

g> n and (7 are the Stieljes transforms of fx n and /i. Moreover, in Lemma 5.5 
in [T^], Schultz proves by using a characterisation theorem of Tillmann that 
l(X) := ir m (L(Al m )) is the Stieljes transform of a distribution A with com- 
pact support in sp(s). Hence, using the inverse Stieljes tranform, Haagerup, 
Thorbjornsen and Schultz get respectively that, for any <p in C£°(R, R), 

- in 

/ipdn n — [ Lpdjj = lim Im [ (p{x)r n {x + iy)dx. (3.12) 
J 71 y^0+ J R 

- in [H] 

ipd/j, n — / (pdfi = lim Im / <p(x)r n (x + iy)dx. (3.13) 

J U 71 y^0+ J R 

Hence, the remainder of the two proofs (in [H] and [Ej) deals with the esti- 
mation of 



lim sup I / cp(x)h(x + iy)dx\ 

where h is an analytic function on C\R which satisfies 

\h(X)\<(\\X\\+KrP k (\Im(X)- l \). (3.14) 

In jH] section 6, Haagerup and Thorbj0rnsen introduce a very clever family 
of functions {I p (X),p > 1} defined by 

1 r+00 

W = 7 TTT / h ( x + 1)^ 1 exp(~t)dt 
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They note that 

h(X) - l[(X) = h(X) 
I p (X)-l' p (X) = I p ^(X), P >2, 
so that for any ip in C£°(R, R) and y > 0, 



(p(x)h(x + iy)dx = / (1 + D) p (p(x)I p (x + iy)dx. 
Jr 

Now, they choose p — k + 1 where k is the degree of the polynomial in the 
right hand side of (j3.14|) (that is p = 8 in and p = 14 in and estimate 
ifc + i(A) for ImX > 0. Using ()3.14|) . it is not difficult to see that 

lim / ^-h(X + z)z k exp(-z)dz = 0. 

r ^+°° J\r,r+ir] kl 



[r,r+ir\ 

Thus, by Cauchy's integral theorem, the authors get 

Ik+x(X) = lim / — h(X + z)z k exp(— z)dz 



1 — >+oo 







J[0,r+ir] k- 

-U(A + (1 + z»(l + i) k+ H k exp(-(l + i)t)dt. 
k\ 



Plugging in ()H.14)) . one gets for any A such that ImX > 0, 

|4+i(A)| < -rr- / (|A| + v / 2t + K) a P fc (|/mA + tr 1 )t fc exp(-t)rft 
- ~kT J (IM + ^t + K^Pkir^expi-^dt 

r+oo 

< / (|A| + v / 2t + i^) a Q(t)exp(-t)rft 
Jo 



fc+1 



where Q{t) = ^jj—P k (t l )t k is a polynomial. 
It follows by dominated convergence 



limsup| / ip(x)h(x+iy)dx\ < / / \(l+D) p ip(x)\(\x\+V2t+K) a Q(t) exp(-t)dtdx < +oo. 
y^o+ Jr Jr Jo 
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Dealing with h(X) = n 2 r n (\) one gets 

f C 
limsupl / (p(x)r n (x + iy)dx\ < — . (3.15) 

Combining (|3.15|) with respectively ()3.12|) and ()3.13|) , one gets respectively 
(fP|) and lOTO) . 



4 The iid case 

We consider a Hermitian matrix X n = [{X n ) k=l °f s i ze 11 f° r which the n 2 
rv ((X n )ii), (v / 2-Re((X n ) ii ) i<i , (\/2i"m((X n )j i ) i<i are independent identically 
distributed with common distribution [ij \fn where fi is a symmetric distri- 
bution with variance 1 onR which satisfies a Poincare inequality (see section 
4.2). We call X n a Wigner matrix with distribution /i. Let Xjp , . . . , Xn be 
r independent copies of X n . We present our main technical tool (see [12]): 

Lemma 4.1 Lei ^ be a real-valued rv such that E(|£| p+2 ) < oo. Let <fi be a 
function from K to C st/c/i £/ict£ i/ie first p + 1 derivatives are continuous and 
bounded. Then, 

E(e0(O) = E^fE(^)(O) + e (4.1) 

(X. 

a=0 

where n a are the cumulants of |e| < Csup t |0( p+1 )(i)|E(|£| p+2 ) ; C depends 
on p only. 

In the following, we shall apply this identity with a function 0(£) given by 
the Stieljes transform of a random matrix. It follows from the Lemma 12.21 
and ()2.12|) above that the conditions of Lemma 14. II fbounded derivatives) are 
fulfilled. 



4.1 The master equation 

Note that since fi satisfies a Poincare inequality, we have j \x\ g dfi(x) < +oo 
for any q in N (see Corollary 3.2 and Proposition 1.10 in ^2]). Note also 
that, since fi is symmetric, any odd cumulant of [i vanishes. 
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Theorem 4.1 With the previous notations, 



E 



J2 a * H nWaiH n (\) + (a - X)H n (X) + 1 

_i=X 

where ||e„|| < P6 ^ Jm ^ A ^ — — and -R n (A) denotes the quantity 



+ -R n (X) + e n = (4.2) 

n 



— E 



^2 ~ ® ln ~ S ^kk a p(^ ® !n - 5 'n)«7 la p(^ ® !» - SrJu a P (A <g> l n 



n 

,p=\ k,l=l 



where k 4 is the fourth cumulant of the distribution \l. Note that 

||^(A)||<P 4 (||/m(A)- 1 ||). (4.3) 



Proof: We shall apply formula (|4.1|) to the .M m (C)-valued function 0(£) 
- Sn)^ 1 for 1 < i,. 



(A (g> l n — S^- 1 for 1 < i,j < n and £ is one of the variable {Xn )u- 



y/2Re((X^') M ), y/2Im{{X^) M ) for 1 < A; < Z < n and p < r. 
We notice that 

90 



dRe((XX 



w - = 4>' x (P)-ekk, l<k<n 



Ikk) 



— — 77 = <t>'( P )-eui, 1 < k < I < n 

dV2Re((X^--^ X " P) 



90 



dV2Im{{X^-^ v 



•'w-fki, l<k<l<n 



kl 



V2Re((X^) kl ) = Tr n (X«e w ), 1 < k < I < n 
V2Im((X^) M ) = Tr n (X^f kl ), \<k<l<n 
(X^) kk = Tr n (X^e kk ), 1 < k < n. 



Let 1 < p < r, 1 < < ? < nbe fixed. For simplicity, we write 0', 0", 0'" for 
the first derivatives of with respect to y/2Re((Xn )ki)- Then, according to 



0' = [(A <S> l n - 5 , n ) _1 a p ® e fe /(A <g> l n — ^n) -1 ^ 
0" = 2 [(A®l n -5 n )^ 1 a p ®e fc/ (A® l n - S^Op ® e kl (\ ® l n -S n y 1 ].. 

0'" = 6 [(A ® l n - S n )~ l a p ® e w (A ® l n - S n ) _1 a p ® e w (A ® l n - S^ 1 
a P <£> e H (A ® l n - S^) -1 ] 
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Writing (J4.1j) in this setting gives 

E[Tr n (X^e kl )(X ® l n - S n )^] = Ie[0'] + + 0(rT 3 ) (4.4) 

where the 0(n~ 3 ) means the norm of this term is smaller than c ^ ap ^ — — . 
Multiplying by n gives the equation, denoted by A^{p): 

nE[Tr n (X^e kl )(X ® l n - S^] = E[<f/] + ^W'} + 0{ n - 2 ) (4.5) 



with the analogous equations with f pq (denoted by B^{p)) and e pp . 

Recall how we can obtain the master equation in the gaussian case (GUE 
case) from (J4.5j) which reads in this case: 

nE[Tr n (X( p )e H )(A®l„-5 ri )^ 1 ] = E [(A ® l n - S n )~ l a p ® e kl (X ® l n - S^' 1 ].. . 

(4.6) 

By a linear combination with the analogous equation with f k i, we have: 
nE[Ti n (X^E kl )(X®l n -S n )^] = E [(A ® 1„ - S^)" 1 ^ ® S«(A <g> l n - 5 n )- x ] .. 
for all 1 < fc, Z < n. 

Now, take in the above formula i = k, j = I and consider -\ ^2 k n we then 
obtain: 

^EpWj^A^ln-S^/] = i [(A ® l n - 5„)^a p (A ® l n - S^ 1 ] 

(4.7) 

that is 

E[id m ® tr n ((l m ® X^)(A ® l n - Sn)- 1 )] = E [tf n (A)a p tf n (A)] . (4.8) 
Now, from the above equation, 



E 



^2a p H n (\)a p H n (X) 
.p=i 



E[i^m ® tr n ((ap ® l n )(l m ® X^)(A ® l n - 



p=i 



£ E[zrf m ® tr n ((a p ® lW)(A ® l n - 5J- 1 )] 
p=i 

E[id m ® tr n (S n - a ® l n )(A ® l n - S^) -1 )] 
_l m + (A - a )E[#„(A)] 
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implying the master formula in the GUE case: 



E 



^a p H n (\)a p H n (\) + (a - X)H n (X) + l r 
.p=i 



0. 



Keeping in mind these computations, we now study the terms coming from 
third derivatives. 

We thus consider A{p) = ^J2k i^fiip) ( res P- B{p)) and study all the con- 
tributions of the different terms. 

Study of the third derivative 

Writing as before the terms appearing in 0'", we can see that all the terms 
except one contains at least two Gki and then, according to Lemma 12.21 
these terms will give a contribution in 0(n~ 2 ) in A(p). The only term to be 
considered is: 

\ E((A <g> l n - S n )ll<X p {\ ®ln~ S^OpiX <g> 1„ - S n )^a p (X <g> 1„ - S^ 1 ). 



n- 



Now, using the same linear combination giving ()4.7|) in the GUE case, we 
obtain that the corrective term of order 1/n appearing in the iid case is: 



1 r K4„ 

-{— E 

n l 2 



1 n 

^ J2 (A ® In - S„)^a p (A ® 1„ - S n )^a p (A ® l n - 



k,l=l 



a p (A ® l n - 5' n ) H 1 ]}. 
The proof of the Theorem is complete. □ 



4.2 Variance estimate 

We assume that \x satisfies a Poincare inequality: there exists a positive 
constant C such that for any C 1 function / : R — > C such that / and /' are 
in L 2 (/i), 

V(/) <cj \ffdpL, 

with V(/) = E(|/ — E(/)| 2 ). We refer the reader to [3] for a characterization 
of the measures on R which satisfy a Poincare inequality (see also pQ). For 
example, fi(dx) = exp(—\x\ a )dx with a > 1 satisfies the Poincare inequality. 
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For any matrices Ax, . . . , A r , define \\(Ax, ■ ■ ■ , A r ) ||g := Y^j=x ll^j'lli- Let 
^ : (M n (C) sa ) r — > W n2 be the canonical isomorphism introduced in Remark 
3.4 in 0. 



Lemma 4.2 For any function f : W n2 — > C be a C 1 -function such that f 
and the gradient V(/) are both polynomially bounded, 

V [/ o . . . , Xf))] < ^E{|| V [/ o . . . , IW)] (4.9) 

Proof := \ij-\fn satisfies the Poincare inequality 

\g- J gd^ n) \W n) <^J \g'W n) . 
f)4.9p readily follows by the tensorisation property of the Poincare inequality. 

4.3 Master inequality 

We follow the lines of the proof of Theorem 4.5 in $,]. Using the master 
equality (|4.2p . we easily get 

|| £[=i a i G n (A)a i G ri (A) + (a - A)G„(A) + l m + ^n(A) + e n \\ 

r 

<\\Y, a *W E {\\Hn(\) -HHn(\))\\ 2 }- 
i=l 

Thanks to (|4.9|) . the following of the proof of Theorem 4.5 in |Hj still holds 
and we similarly get 

E{||F n (A) - E{H n {\))\\ 2 } < -5-H 5> 2 IIII(MA))-T 
and therefore 

|| Va,G n (A)a l G n (A) + (a - A)C n (A) + l m + --R n (A) || < MIKM^) !J!1 , 

(4.10) 
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4.4 Estimation of \\G n — G\\ 

In the Gaussian case, Haagerup and Thorbj0rnsen in [§] and Schultz in 
prove that G n (X) is invertible for any A such that ImX is positive definite. 
In our more general case, we are going to use the master inequality (j4.10j) in 
order to prove that, for any A in some subset of O, G n (X) is invertible and 
to get an upper bound of ||G n (A) _1 ||. 
Set 

r 

B n (\) = ^2a i G n (X)a i + (o -A). 
i=i 

Now, from the master inequality (|4.1U|) and (|4.3j) . we get 

||^a l G n (A)a J G„(A) + (a - X)G n (X) + l m || < MK^B> . ( 4 .n) 

i=i 

that is 

i, B „ (A)G „ (A) + M < mum-id 



n 



Hence, for any A such that Pe ^ Im ^ — HI < I ( B n (X)G n (X) is invertible with 

||[ J Bn(A)G n (A)]- 1 || <2. 
Thus, for such a A, G n (A) is also obviously invertible with 

r 

HG^A)- 1 !! < 2||5 n (A)|| < 2(||a„|| + ||A|| + £ ||a i || 2 ||(Jm(A))- 1 ||) (4.12) 

i=i 

Now, from the inequality (|4.11|) and using ()4.12|) . we get readily that for any 
A in O such that MMTl < i 

n 2 ' 

|| <kG n (\)oi + (a -A)+G n (A)- 1 || < P 6(||(MA)) 2(||ao||+||A||+^ |k|| 2 ||(/m(A))- 

i=l i=l 

(4.13) 

Define 

o' n = {xeo, MKMA)nD 2(||ao|| + ||A|| + £ || aif ik^))-!,,) < 1 



- u " 11 11 ^" *" - ^ " 2||(/m(A))-i| 
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Since t i— > ^— ^ — -2(||ao|| + £ + 21=1 ll a «ll 2 ^ 1 is a continuous strictly 
decreasing function from ]0, +oo[ onto ]0, +00 [, one can prove that O n is an 
open connected subset of M m (C) by following the proof of (a) Proposition 
5.6 in jUJ. Note that, using the inequality 



< 



||(/m(A))-i| 
one immediately gets that for any A in 0' n , 

r 

||(/m(A)) -1 ||(||ao|| + ||A|| + £ ||ai|| 2 ||(Jm(A)) -1 ||) > 1 

1=1 

and thus that 

PedKMA))- 1 !) < 1 
n ~ 4' 

Consequently, for any A in 0^, G n (X) is invertible and (|4.13|) holds. Defining 
for any A in O n , 

r 

An (A) = a + J] a * G «( A ) a i + Cn(A)" 1 

i=l 

and sticking to the proof of P described in section II. 2.1, we get that, for 
any A in O, 

\\G n (X) - G(X)\\ < 4 p 6(||(MA))- 1 ||) + + A || i^ikj^j-iidikj^))- 

n 

< ( || A || +Jf) ^(IKMA))-||) 

n 

Note that, in the following we will use (|4.13|) in the simplest form: 
W\eO' n , || Va i G n (A)a i + (a -A)+G n (A)- 1 || < (\\X\\+K) MK^lUl. 

4.5 Convergence of R n (X) 

Let Xj, i = 1, . . . , r be self-adjoint operators in a C* probability space (£>, r). 
We assume that the Xi are free and identically semi-circular distributed with 
mean and variance 1. Then, G satisfies (|3.1|) . 
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Proposition 4.1 Let a be a matrix in M m (C). Then, 



E[i 5>®1„-S„)j£ a (A®l n -^)^] = GaG+O ( Pl ° (ll(MA)) " 1|l) (||A|| + #)) 
fc=i ^ ' 

(4.15) 

Proof: We start from the resolvent identity: 

n 

A(A <g> 1„ - S^)^ = l m + y^(<S'n)fc;(A <8> l n - S'n)^ 1 - 

i=l 

We write G^(X) = (A®l n -S' n ) _1 and L>i n) (A) = ± ELi ^(A)^ a G^{X) kk . 
From the above identity, 

1 n 

A-VG?W(A) fcfc oG^(A) fcfc 

n n 

a-y^G^(X) kk + - y2(S n ) kl G^(X) lk a G (n \X) kk 
fc=i fc,i=i 

1 r n 

atf n (A) +a 0J Di n) (A) + -^a p ^(lW) H G (n) (A) tt a G^(A) 



n 

k=l 



kk 

n — ' * — ' 

p=l k,l=l 



We take the expectation and we use the integration by part formula ()4.1|) for 
the last term: 

E((X^) k MX'^ .Y<;»)] = EfIY„(A'«'»B It )*(.f . . . , A'")] 

= iE[ W ,,,,,^,),^ 1+0 ( ft(ll(W-||) 
n \ rr 

with $(X«, . . . ,xl r) ) = GW(A) tt a GM(A) fcfc . Then, 

E[(XW) JH G?( n )(A), fc aG?W(A) fcfc ] = -E[G(")(A) ii a p (?W(A) tt aGW(A)^] 

n 

+^E[G^(X) lk aG^(X) kl a p G^(X)i k ] 



+0 



P 4 (||(Jm(A))- 1 ||) 



n 2 
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Thus, we obtain from the resolvent identity, 
(X-a )E(D^(X))=aE(H n (X)) + 

r n 1 

Y^a p — E[G (n) {X) u a p G {n) (X) kk aG {n) (X) kk + 



p=l fc,Z=l 



GW(A) ft aGW(A), ap G(")(A),]+o( f4(ll(Jm w (A)) " 1|l) ] . 



^From Lemma l2. 2| 

1 £ E[GW(A) ifc aGW(A) H a p GW(A) lfc ] = O ^MM^W) 



fc,/=i 
thus, 

r 

(A-ao)EpW(A)) = aE(F„(A))+^a p E[if n (A)a p DW(A)]+0 

P =i 

From the estimate of the variance of ff„, we have: 



P 4 (||(/m(A))-l) 



n. 



E[fl n (A)a p £)W(A)] = E[H n (X)]a p E[D^(X)) + O 



p^IKMA))- 1 !!) 



Then, using also the estimation of ||G n (A) — G(A)|| we get 
(\-a Q -J2a p G(\)a p )E(D^\\)) = aG + O f ^M^W) m + K) \ 

Using (|3.1|) and ()2.9j) we finally get 

E(^(A)) = GaG + ^(IKM*))" 1 !!) (||A|| + K) \ . 

□ 

^From the above proposition, we obtain: 

Proposition 4.2 R n (X) defined in Theorem \4-l\ converges as n tends to in- 
finity to 

r 

R(\) = ^-^2a p G(X)a p G(X)a p G(X)a p G(X). 
P =i 
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More precisely, 



\R n {\)-R{\)\\<{\\\\\+K) 



2- P 20(||(-MA))~ 1 ||) 



n 



(4.16) 



Proof: It's enough to prove the convergence of each coordinate of the mxm 
matrix i? n (A). This will actually follow from the convergence of terms of the 
form: 



E 



a,/3 



n 



kk 7,<5 \ n 



- 1 



k=l 



k=l 



(4.17) 



for a et b elements of the canonical basis in A4 m (C). Since, applying Lemma 
14.21 we have 

E(||I?W(A) - EpW(A))|| 2 ) < P6(ll(Jm(A)rl|l) ; 



the above quantity (|4.17jl is of the same order as: 



E 



«,/3 n 



-i r«( n ) / ,r<( n ) 

2_^ U kk aLr kk 



k=l 



E 



7,<5 ( n 1 G kk bG kk 
k=l 



According to Proposition ^. 1| this last quantity converges towards a ^{GaG) 1: §(GbG). 
Thus, the convergence of R n to R follows with the estimation ()4.1fij) . □ 

We define 

L(A) = {id m ® r)[(A ® 1 B - s)- 1 (i?(A)G(A)- 1 ® 1 B )(A ® 1 B - s)" 1 ]. 

4.6 Estimation of ||G(A) - G n (A) + ±L(A)|| 
Following 4.24 in (THj, one gets for any A in 0' n , 

||G(A)-G n (A) + ±L(A)|| 

< || (A ® 1 B - s)- 1 !! ||(A n (A) ® 1 B - s)- 1 !! ||A„(A) - A + --R(A)G(A) _1 || 

n 

+-||(A ® Is - s)- 1 ||||i?(A)G'(A)- 1 ||||(A ® l fl - s)- 1 - (An (A) ® 1 B - 1 1| 
n 

< ||(/m(A))- 1 ||||(/m(A n (A)))- 1 ||||A n (A)-A + - J R(A)G(A)- 1 || 



n 



c, 



+-||(/m(A))- i |n|(/m(A n (A)))- 1 ||||A„(A)-A||(||A||+K)^ 
n 
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where we made use of the estimates ()2.7|) . (|2.8|) . (|2.6|) and the upper bound 

II^hcikma))- 1 !! 4 . 

Now, for any A in O n , 

||(/m(A n (A)))- 1 ||<2||(/m(A))- 1 || 

and 



\KW->~\\< m{Im{X)) ^\\W\ + K). 

n 



Thus, 



\\G{\)-G n {\) + -L{\)\\ < 2||(/m(A))- 1 || 2 ||A n (A)-A + - J R(A)G(A)- 1 | 
n n 



^ 14 (ll(MA))-'||) (||A||+g) 3 



Now, for any A in O n , 



||A n (A)-A + -i?(A)G(A)- 1 || < ||A n (A)-A + - J R n (A)G„(A)- 1 || 
n n 

+-|| J R n (A)(G n (A)- 1 -G(A)- 1 )|| 
n 

+-\\(Rn(X) - R(X))G(X)-')\\ 
n 

< ^(ll(MA))-'ll) 

n 2 

+ ft(ll(MA))-'||) . 

n 

+i|| J R n (A)- J R(A)||(||A||+K) 2 ||(/m(A))- 1 | 

where we used (l4.1Uj) . (|4.12j) . (|4.3|) and (|2.8|) . Moreover, one easily gets 

||G n (A)- 1 -G(A)- 1 || = ||G n (A)- 1 (G(A)-G n (A))G'(A)- 1 || 

< ||G n (A)- 1 ||||G(A)-G n (A))||||G(A)- 1 ||. 

Consequently, using the estimate (j4.14|) of ||G n (A) — G(A)|| together with 
gUP) and <fgfg|i - we get 

|| G „(A)-- G (A)-||<(||A|| + A-r P "< ll < /m(A) >" ll) . 

n 

23 



We conclude that 

P ir (\\(Im(X))- 1 \\) 

2 



||G(A)-G„(A) + -L(A)|| < (\\X\\ + K) 

n n z 



+-\\R n (\) - R(\)\\(\\x\\ + miiiMm 



- 1 1 1 3 



n 



Using ()4.16jh we can conclude that, for any A in O n , 

||G(A) - G„(A) + ii(A)|| < (||A|| + ^4 ^(ll(/m(A))-'||) 
Now, for A in 0\0'„, 

1 < 4 Pe(ll(Jm n (A)) " 1|l) (|| a „|| + ||A|| + || g a ?||||(/m(A))-'||)||(Jm(A))-'||} 

< (||A|| + Jf) P 8 (ll(MA))-||) 

n 

We get 

||G(A)-G n (A) + -L(A)|| < ||G?(A)-G?„(A)||+-||L(A)|| 
n n 



n 

x r .. /v 



A |l + g) MKM^ + i ||(Jm(A)rl||7(||A||+ ^ 



< m+K) 



n n 
P 17 (\\(Im(X))-i\\) 



n 2 



Thus, one can easily see that one can choose K and P23 such that for any A 
in 0, 

||G(A) - G n (X) + ±L(A)|| < (||A|| + ^ ^(IKMA))- 1 !!) ^ 

Note that, since under our hypothesises, S n and — S n are identically dis- 
tributed, the arguments of to prove her theorem 4.5 still hold. Thus, 
f!4.18|) is also valid for any A such that ImX is negative definite. 
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4.7 Spectrum of S n 

• From step 2 to step 3 

Sticking to the proof of Lemma 5.5 of fT5| . we get that, 

/(A) :=tr m (L(Xl m ),X G C\R, 

is the Stieljes transform of a distribution A with compact support in sp(s). 
Hence, the proof described in section 3.2 still holds (with a = 4 and k = 23); 
thus we can state that for any smooth function (p with compact support 

E[(tr m ®tr„)(</?(£ n ))] = (tv m ®T)(<p(s)) + -A(<p) + 0(— ). (4.19) 

Moreover, following the proof of Lemma 5.6 in [To] , one can show that A(l) = 
and deduce that, for (p smooth, constant outside a compact set and such 
that supp(ip) fl sp(s) = 0, 

E[(tT m ®tT n )(<p(S n ))] =0(i). 

• Step 4 

The proof of step 4 is exactly the same as in jH] so that we have proved that, 
for any e > and almost surely 

Spect(S n ) C Spect(s) + (— e, e) 

when n goes to infinity. Note that this result implies that 

s«p n ||X^|| < +oo a.e. 

4.8 The main theorem 
4.8.1 First inequality 

By the same arguments of in Proposition 7.3, we can deduce the following 
inequality from the above inclusion of the spectrum of S n . 

Proposition 4.3 Almost everywhere, for all polynomials p in r non com- 
muting variables, 

limsup|b(X«,...,xW|| < \\p( Xl ,..., Xr )\\. 

n— »+oo 
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4.8.2 Second inequality 

Proposition 4.4 Almost everywhere, for all polynomials p in r non com- 
muting variables, 

liminf|b(XW,...,4 r )|| > \\p( Xl ,...,x r )\\. 

n— >+oo 

^From [S], Proof of Lemma 7.2, it is clear that this proposition follows from 
the almost sure asymptotic freeness of the Xn together with the property 
that 

a.esup ll-X^H < +oo. 



it 



The proof of the first point follows the proof of Theorem 6.2 in [To]; never- 
theless, we modify the proof of Lemma 6.5 in 15 to get the analogue in our 
context without needing such a result as Lemma 6.4 in |15j . 

Lemma 4.3 Let d be in W , i ly . . . ,i d be in {1, ... , r} and n be in N*. Define 
f : M n (C) r ^Cby 

f(v t , ...,v r ) =ti n (v il ...v id ). 
Then, there is a constant C > such that 

C 



V/(lW...,lW)< 



n 2 



r 
sa 



Proof: Applying Poincare Inequality (|4.°-j) . we get 

v/(x B w, . . . , XW) < £e{||v/(xW, . . . , X«) II*}. 

Now, let v — (t>i, . . . , v r ) be in M n (C) r and w = (wx, . . . , w r ) be in M n (C) 
with || w || e = 1. By the Cauchy Schwartz inequality 

1^- /(« + ttu)| = -\Trn(w h v i2 v i3 ...v id ) + Tr n (v h w i2 v h ...v id ) + ... 
att=o n 

... + Ti n (v h v i2 . ..v id ^w 

< - {IKilWka^a • ■■ v idh + ■■■ + IKJWki^ ■ ■ •Ui.j.Jla} 

< - {IK 2 "Ui 3 • ••fiJh + H^ii^s • ■ -^Jh + • • • + ll^i^ia • • -Vid-ih] 
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1-1 2 

2 



Thus, 

+ . . . + {Tr n (X^ . ..x^X%+-4 . ..Xt ] )} 

< - 2 {Tr n (X^X^ . . . XM X M . . . Xg*) 

+ ... + Tr n (X™ . . . Xt^X^ . . . X™)} 

for some constant C depending on d, and we get that 

E(||V/(X«,...,XW)|P < ^{E(tr n (Xt ) Xt ) ...X^X^...Xt ) )) 

+ ... + E (tr„(X^ . . . X^X^ . . . XW)) } . 

Each term inside the brackets of the left hand side is uniformly bounded in n 
since it converges as n tends to infinity according to the result of asymptotic 
freeness in mean of Dykema in |]] . The result follows. □ 

Lemma 14.31 yields the almost sure asymptotic freeness of the Xn using 
the Borel Cantelli lemma. 



In conclusion, 

Theorem 4.2 Let Xn \ • • • , X^p be independent Wigner matrices associated 
to a symmetric distribution fi which satisfies a Poincare inequality. Let 
(x\, . . . , x r ) be a semicircular system. Then, almost everywhere, for all poly- 
nomials p in r non commuting variables 

lim \\p(Xi 1 \...,X^)\\ = \\p(x 1 ,...,x r )\\. 

n—>+oo 



5 The Wishart case 

We consider a n x n Hermitian matrix Y, distributed as a Wishart matrix of 
parameter p(n) > n and variance - that is with density w.r.t the Lebesgue 
measure dM on A4 sa (C): 

C n , p l ( m>o) (det (M) f^- n exp(— ntr(M)). 
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V[ TL ) 

We assume that ► a for some a > 1. More precisely, according to 

71 n— >oo 

Dirichlet theorem Q16J, Lemme 14.1), there exists subsequences p(n) and 
q(n) of integers tending to oo such that: 

\P{ n ) i ^ 1 
— a < 



g(n) g(n) 2 

So, we shall consider a matrix F of size q(n) and parameter p(n). For sim- 
plicity, we shall denote the subsequence q(n) by n and therefore, we will 
assume in this section that: 

a < — . (5.1) 

n rr 

It is well know that the spectral measure of Y converges to the so called 
Marchenko-Pastur distribution \i a [T3] : 

Va{dx) = — l [(VS _ 1)a)(V5+1)a] (x)da;. 

5.1 Differentiation formula for the Wishart ensemble 

Lemma 5.1 Let $ a C 1 function on A4 sa (C) tyit/i $(0) = ; then: 

E[$\Y).H]-nE[$(Y)Ti n (H)] + (p(n)-n)E[<3?(Y)Tr n (Y- 1 H)] =0 (5.2) 
for all hermitian matrix H , or by linearity for H = E^, 1 < j,k < n. 
Proof: Since the Lebesgue measure is invariant by translation, 

E[$(y)] = y $(M+eX)exp(-nTr n (M+eX))(det(M+eX)) p W~"l (M+eX > 0) dM. 

Now, by differentiation with respect to e and taking e = 0, we obtain ()5.2j) 
using <9(det M) s = s(det M) s M~ l . 

5.2 The master equation 

Let (Xn i ■ ■ ■ , Xn^)i=i,... r be r independent copies of the random matrix Y. 
We shall apply (|5.2j) with 

$(X«) = [(l m ® X«)(A ® 1„ - Sn)' 1 ] jk G A4 r 
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and H = Ej k . Then, 

&{X®).E jk = (A®l n -5 n ) fefc 1 +[(l m ®X«)(A®l n -S , n )- 1 ] J ^[(A®l n -S„)- 1 ] fefc 
and 

$(Xg>) Tr n ((XV)-% k ) = (Jf (0)^? [(l m ® X«)(A ® l n - ^J- 1 ] jk . 
The sum over j of the terms in the above equation gives: 

(A ® l n - Sn)^ 1 . 

Now, if we sum the identities obtained by (|5.2|) over j, /c, and dividing by n 2 , 
we obtain: 

E[id®tr n ((A® In-Sn)' 1 )] (5.3) 
+E[id ® tr n ((l m ® X^)(A ® l n - S^ 1 ) a, id ® tr„((A ® l n - S^ 1 )] 
-E[id ® tr n ((l m ® X«)(A ® l n - S 1 ,,)- 1 )] 

+ (W _ i) E [^ g, trn (( A g i n _ Sn)- 1 )] = 
which can be written as: 

E [id ® tr n ((l m ® X^)(A ® l n - S n )~ 1 )(l m - a, id ® tr n ((A ® l n - S,,)" 1 )] 

= ^E[*d <g> tr n ((A ® l n - S^)" 1 )]. (5.4) 

Proposition 5.1 1. For A 6 O, 

|E [id ® tr n ((l m ® X®)(\ ® 1„ - S n )- 1 )(l m - a z id ® tr n ((A ® l n - S^)" 1 )] - 
E [id ® tr n ((l m ® )(A ® l n - Sn)" 1 )] E [(l m - a, zd ® tr n ((A ® l n - S n )' 

< P 4 (||(/m(A))- 1 ||) 
— n z 

£ For AGO, 

HG^Ani^diAii + KflKMA))- 1 !! 
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3. If ai is invertible and A G O, then (l m — aiG n {\)) is invertible and 

IKU-cnGnCA))- 1 !!^ Ha^lKHAii + ^iKMA))- 1 !!. 

J/||(/m(A)) -1 || < 2pnT' ^ en (l m ~ a iCn(A)) invertible and 

\\(l m -a l G n (X))- l \\<2. (5.5) 

Sketch of Proof: 

1 . The variance estimate follows from the Gaussian Poincare inequality since 
we can write Y = ^X*X for a rectangular Gaussian matrix X. We proceed 
as in j3 Section 4]. We need some estimate on the maximal eigenvalue of 
Y, i.e. E[A max ] and E[A^ a J are bounded, independently of n. This can be 
proved, as in Lemma 5.1 of 9], using previous results in [S]. 

2. The proof is the same as Proposition 5.2 in 

3. If ai is invertible, 

(1 - aiG n (X)) = ai{a{ 1 - G n (X)) := a x T. 

Now, the matrix T satisfies, Im(T) = —Im(G n (X)) and thus is positive 
definite (see the proof of Proposition 5.2 in [§]). Its inverse T -1 satisfies: 

HT- 1 !! < HMT)- 1 !! = H/m^A))- 1 !! < (||A|| + ^) 2 ||(Jm(A))- 1 ||. 

(15.5)1 follows from the majoration ||G n (A)|| < ||(Im(A)) -1 ||. □ 

From (|5.4|1 and Proposition 15. 11 we obtain: 

r 

a{E[id ® tr n ((l m ® Y®)(\ ® l n - S' n )- 1 )]- 



p(n) 



n 



^a^E[iJ„,(A)](l m -a i E[^„(A)])- 1 



<(l|A||+#) 



PsdKMA))- 1 !!) 



n 2 



The first line of the above equation equals: 

E[id ® tr n (S n - a <g> 1„)(A ® l n - S^ 1 )} = -l m + (A - a )E[fl n (A)]. 
We have thus obtain the master inequality: 

||-a G n (A)+AG n (A)-^Va / G n (A)(l m -a / G n (A))- 1 -l m || < (||A||+#) 2jFM||( //,,(A)) 



n *■ — ' n 2 

i=i 



(5.6) 
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or since the matrices (l m — aiG n (X)) 1 and aiG n (X) commute, 

||-a G n (A)+AG n (A)-^V(l m -a / G n (A))- 1 a z G n (A)-l m || < (\\X\\+K) 2mil '" {X)) 



n *■ — ' ' n 2 



1=1 



(5.7) 



5.3 Estimation of \\G n (X) - G{\)\\ 



Let X{,i < r be a free family of self adjoint variables in a C*-probability 
space (£>, r), with Marchenko-Pastur distribution fi a , with parameter a. 

r 

s = a Q <g> l s + 2J a « ® 

i=l 

Using the known expression of the R transform of the distribution of Xi (see 
0, [El Example 3.3.5] 2 ): 

R x {z) = a(l - z) -1 ,* e C\R, 
we can show the following 

Lemma 5.2 G satisfies the following equation: for X e O, 

r 

a + a ^2(l m ~ OiGiX))- 1 ^ + G(A)- 1 = A. (5.8) 
i=i 

Sketch of Proof: From the definition of the R transformation with amal- 
gation over M m (C), we can show that: 

RamW = Rx(aX)a, a 6 M m (C) sa , A G O 
and then, by freeness asumption 

r 
i=l 

()5.8|) follows, using the relation between R and G. □. 



2 We warm the reader that the R transform defined in this book differs by a factor z 
from the Voiculescu R transform we used here 
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Theorem 5.1 For any A G O, 



P 8 (||(/m(A))- 1 |i) 



||G(A) - G B (A)|| < (||A|| + Kf — X " (5.9) 



n 



Proof: We can proceed as in the proof of Theorem 5.7 in [§]. We just 
mention the different steps: 

Step 1: Define A n (A) = a + G„(A) _1 + a ^(l m — aiG n (X))~ 1 ai. From the 
master inequality (J5.7|) . Proposition 15 .11 and (|5.1j) . we can show that: 



.PedKMA))- 1 !!) 

Then, for A 6 Ol of the form 



|A-A n (A)||<(||A|| + ^ 

n z 



M\\{Im{\))^\\) ^ 1 



0' n = {\eO, (||A|| + Kf HHlv v " »" < 



n 2 2||(/m(A))- 1 | 
we have Jm(A n (A)) > 2||(jm(A))- 1 || an< ^ m particular A„(A) G O. 

Step 2: For A G C^, we can consider G(A n (X)) and we have, from the 
identity (Q: 

a +G n (A)" 1 +a ^(l m -a z G n (A))- 1 a / = a +G(A n (X))- l +a ^(l m -a i G(A n (A)))" 1 a i . 

(5.10) 



Lemma 5.3 (see Propostion 5.6) For A G 0' n , 

G(A n {X)) = G n {\). (5.11) 

Proof: As in j^j, it's enough to prove (|5.11|> for A G O", a non empty subset 
of the connected subset 0' n . Put x = G n (X) and y = G(A n (A)), then, from 
ffOTIjl . 

a + a;" 1 + a 2j(l m - a/a;) -1 ^ = a + y _1 + a 2j(l m - cny)~ l ai 
i i 

so that 

?/ + a z(l m - aix)" 1 ^ y = x + a ^ x(l m - aiy)~ x ai y, 
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Thus, 

y-x = ay^x[(l m - aiyY 1 ~ O-m - aix)' 1 }^ y 
i 

= x 0-m - o-ivY 1 ^ - a i x ) ~ (1 - aiy)]{ l m - aja:) -1 ^ y 

i 

= a ^2 x ( l rn - aiy)~ x ai{y - x)(l m - aix)" 1 ^ y 
i 

In particular, we have, 

\\V ~ X \\ < ^IMIIMI ^2 ll(lm - aiy) _1 ||||(lm - a/x) _1 ||||a;|| 2 ^ - x\\ 

(5.12) 

Now, we have 

NMI^A)!! < ||(/m(A)ni 

and 

||y|| = ||G(An(A))|| < (||(/m(A n (A))) _1 || < 2||(/m(A))- 1 || 
for A G Q' n (see Step 1). 

Moreover, from Proposition 15. 1| for ||(Jm(A)) _1 || small enough, 
IK^-a,*)- 1 !! <2; IK^-a^)- 1 !! < 2. 

Set 

0£ = {A G O^; allGf^HIIGfCA^CA))!! IK^-ajGCA^CA)))- 1 !!!!^^^^))- 1 !!!^,!! 3 < 1}, 

z 

Then, from (pTT2j) . for A G 0£, G(A n (A)) = G n {\). Now, it is easy to see, 
from the above estimates, that A = itl m G O'^ for t large enough, so O'^ is a 
non empty set. □ 

Step 3: The estimation of G(A n (A)) —(7(A) is obtained as in Subsection 3.1 
(considering the two cases A G 0' n and A G C?\C^J. □ 

5.4 The spectrum of »S n 

^From Theorem 15. II and the proof described in Section 3.2 (see also Section 
6 in 9j), we can prove that for (p smooth, constant outside a compact set 
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and such that supp(<p) fl sp(s) = 

E[(tT m ®tv n )(<p(S n ))} = O(-). 

from which we deduce that, for any e > and almost surely 

Spect(S n ) C Spect(s) + (—e,e) 
when n goes to infinity. 

5.5 The main theorem 

We can now prove: 

Theorem 5.2 There exists a set N of probability such that for all non 
commutative polynomial p in r variables, and all u e Q\N , 

lim | |p(X«0"), X ( n r) (uj))\ | = | \p( Xl , ...x r )\\. (5.13) 

n — ► oo 

Proof: The inequality 

lim sup | \p{X^ ,...,X^)\\ < | \p( Xl , ...x r )\\ a.s. (5. 14) 

n ► oo 

follows from the above inclusion of the spectrum of S n and the arguments 
developed in [Hj, Section 7. The reverse inequality 

liminf ||p(X«,...,X«)|| > ||K^i,...x r )|| a.s. (5.15) 

n ► oo 

follows, as in Lemma 7.2 in pj, from the a.s. asymptotic freeness of the 
(Xn )i=x t ... } r and sup n oo a.s.. The first point was proved by Hiai 

and Petz (see ^UJ , ^1]) and the second point follows from ()5.14|) . □ 

Remark: If we only assume the convergence of to a with — a\ < -, 

° n 1 n 1 — n ' 

then an extra term appears in the estimation of G — G n at order n~ 2 , namely: 
\\G(X)-G n (X) + (a-^-)L(X)\\<^- 

with 

L(X) = (id m ® r)[(A ® 1 - s) _1 (fl(A) <g> 1 B )(A ® 1 B - s)" 1 ] 
and .R(A) = J];! 1 ™ - a;G(A)) -1 a/,. 

As in Schultz [T3] and in the iid case (see Section 4), this term gives rise to a 
distribution with compact support in sp(s) and the conclusion remains true. 
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